Semileptonic A& >c to Nucleon Transitions in Full QCD at Light Cone 



> 
O 



K. Azizi 1 ^, M. Bayar 2 **, Y. Sarac 3 '*, H. Sundu 4 <* 
t Physics Division, Faculty of Arts and Sciences, Dogus, University, Acibadem-Kadikoy, 

34722 Istanbul, Turkey 
* Department of Physics, Kocaeli University, 41380 Izmit, Turkey 
* Electrical and Electronics Engineering Department, Atilim University, 06836 Ankara, Turkey 

1 e-mail:kazizi@dogus. edu.tr 
2 melahat. bayar@kocaeli.edu.tr 

3 ysoymak@atilim. edu.tr 
4 hayriye. sundu@kocaeli.edu.tr 

On ! 

The tree level semileptonic At — > plv and A c — > nlv transitions are investigated using the light 
cone QCD sum rules approach in full theory. The spin-1/2, Aq baryon with Q — b or c, is 
£\j ' considered by the most general form of its interpolating current. The time ordering product of 

the initial and transition currents is expanded in terms of the nucleon distribution amplitudes with 
different twists. Considering two sets of independent input parameters entering to the nucleon 
wave functions, namely, QCD sum rules and Lattice QCD parameters, the related form factors and 
I , their heavy quark effective theory limits are calculated and compared with the existing predictions 

of other approaches. It is shown that our results satisfy the heavy quark symmetry relations for 
I ' lattice input parameters and b case exactly and the maximum violation is for charm case and QCD 

sum rules input parameters. The obtained form factors are used to compute the transition rates 
i— i . both in full theory and heavy quark effective theory. A comparison of the results on decay rate of 

p£ ' At — > plv with those predicted by other phenomenological methods or the same method in heavy 

Oh' quark effective theory with different interpolating current and distribution amplitudes of the At is 

^ , also presented. 

<D ' 

f-H 1 PACS numbers: 11.55.Hx, 13.30.-a, 14.20.Mr, 12.39.Hg 

<N ■ I. INTRODUCTION 

> ' 

oo : 

, Motivated by the recent experimental progresses on the spectroscopy of the heavy baryons containing heavy b or c 
r**- ■ quark [J [H, H, EL IE IE 0> Hi > theoretical studies on these baryons gain pace. Because of the heavy quark, these states 
are expected to be narrow, experimentally, hence their isolation and detection are easy comparing with light systems. 
Theoretically, investigation of the semileptonic decays of the heavy baryons, whose experimental testing may be in 
the future program of the large hadron collider (LHC), have attracted interests beside their mass and electromagnetic 
ON . properties. For instance, the semileptonic Ah — > A c and A c — > A decays have been investigated in three points QCD 
sum rules and heavy quark effective theory (HQET) in 0. The Ah — > plv transition has also been studied in the same 
frameworks in [l(| and using SU(3) symmetry and HQET in 11]. Constituent quark model have also been used to 
study the A c — > nlv and Ah — ► plv form factor s 11211 and semileptonic decays of some heavy baryons containing single 
heavy quark in different quark models (l2l . Il3l [Tj] are some other works in this respect. 

In our recent work (l5j . we analyzed the semileptonic decay of Eh, which has different interpolating current and 
structure than Aq, to proton in light cone QCD sum rules. In present study, we calculate the form factors related to 
the semileptonic decays of the Ah — > plv and A c — > nlv also in light cone QCD in full theory and HQET limit. In full 
theory, these transitions are governed by six form factors, but heavy quark effective theory limit reduces them to two. 
The vacuum to nucleon matrix element of the time ordering product is expanded in terms of nucleon distribution 
amplitudes (DA's) near light cone, x 2 ~ 0. The nucleon wave functions contain eight inde pen dent parameters, which 
we consider two sets, namely, calculated using the QCD sum rules [l6j] and lattice QCD jl7l . [l8l Il9| approaches. In 
the calculations, the most general current of Aq generalizing the Ioffe current is used. The obtained form factors 
are used to compute the corresponding transition rates both their numerical values and in terms of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements. Studying such type of transition provides a better understanding of 
the internal structure of Aq, information about the DA's and input parameters as well as determination of the CKM 
matrix elements. Note that, using different interpolating field, the semileptonic decay of bottom case, Ah — ► plv, 
has already been investigated in references 0, [2l[ in the same framework but HQET limit. In [2(|, the nucleon 
distribution amplitudes are used only with QCD sum rules input parameters, while the distribution amplitudes of Ah 
have been utilized to calculate the form factors in [21] . 

The layout of the paper is as follows: in section II, the details of the calculation of the form factors in light cone 
QCD sum rules method are presented where the nucleon distribution amplitudes and the most general form of the 
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interpolating currents for the Aq baryon are used. The heavy quark limit of the form factors and the relations between 
the form factors in this limit is also discussed in this section. Section III comprises numerical analysis of the form 
factors and our predictions for the decay rate obtained in two different ways: first, using the DA's obtained from QCD 
sum rules and second, the DA's calculated in lattice QCD. A comparison of our results on form factors and transition 
rates with the existing predictions of other approaches is also presented in this section. 



II. THEORETICAL FRAMEWORK 

In this section, following fl5j |. we calculate the form factors of the Af, — > p and A c — > n transitions in the framework 
of the light cone QCD sum rules and full theory. At quark level, these decays are governed by the tree level Q — > q 
transition, where Q represents b (c) quark and q stands for u (d) quark for Af,(A c ). The effective Hamiltonian 
responsible for these transitions at the quark level has the form 



^V qQ q^(l - 75 W(l ~ J5>- (1) 

To calculate the amplitude, we need to sandwich the above equation between the initial and final states and compute 
the matrix element (iV| ^-/^ (1 — 7 5 )Q|Aq), which is needed to study the Aq — > Nlv decay. The starting point is to 
consider the following correlation function 

IL>,g) = i ( d 4 xe^(N(p) \ T{ J**(x) J A « (0)} | 0), (2) 



where, J Aq is interpolating currents of A b( - C ) baryon, J^ r = <77 M (1 — J5)Q is transition current and (N(p) | represents 
the nucleon sate, where p denotes the proton (neutron) momentum and q = (p + q) — p is the transferred momentum. 

One further step of the calculation is the saturation of the correlation function by a tower of hadronic states having 
the same quantum numbers as the interpolating currents. The obtained result from this procedure is called the 
phcnomenological or physical side of the correlation function. From the general philosophy of the QCD sum rules 
approach, this correlator is also calculated using the operator product expansion (OPE) in deep Euclidean region. 
This part is called the theoretical or QCD side. Match of these two different representations of the same correlation 
function gives sum rules for form factors. To suppress the contribution of the higher states and continuum, the Borel 
transformation is applied to both sides of the sum rules for physical quantities. 

Let first calculate the phenomenological part. After the insertion of the complete set of the initial hadronic state 
and performing the integral over x, we obtain the physical side as: 



_ (N(p) | J*[{x) | A Q (p + q,s)){A Q (p + q,s) | J A «(0) | 0) 



IL>, q) = ^ ^ 'V 2 + -> ( 3 ) 

where, the ... represents the contribution of the higher states and continuum. The matrix element (Aq(p + q, s) \ 
J A « (0) I 0) in © is given by: 

(A Q (p + q, S ) I J A «(0) \0)=X AQ u AQ {p + q,s), (4) 
where Aa q is residue of Aq baryon. The transition matrix element, (N(p) \ J 1 ^ \ Aq(p + q, s)) can be written as 



(N(p) I J l ;{x) I A Q (p + q)) = N(p) [ 7m /i(Q 2 ) + ™»Xh{Q 2 ) + q^hiQ 2 ) + 7m753i(Q 2 ) + ^^q v g 2 {Q 2 ) 



<Z M 75<? 3 (Q 2 ) 



ua q (p + q), 



(5) 



where Q 2 = —q 2 . The fi and gi are transition form factors in full theory and N(p) and ua q {p + q) are the spinors of 
nucleon and Aq, respectively. Using Eqs. ([3]), (UJ) and and summing over spins of the Aq baryon, i.e., 



u ^q (p + 9; s W Q (p + q> s ) =0+ 4 + to a« i 



(G) 
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we attain the following expression 
Aa 



n M (p,<?) = 



(p + qf 



N(p) [7m/i(0 2 ) + i^uq"f2(Q 2 + q^hiQ 2 ) + 7m753i(Q 2 ) + ^^q v g 2 (Q 2 ) 



4+ d + m AQ ) + 



Using 

Na^q^UAQ = i N[(m N + toa q )7p - ( 2 P + q)^}uA Q , 
in Eq. ([7]), the following final expression for the physical side of the correlation function is obtained: 



(7) 
(8) 



Ha(p, q) 



m A Q - + <?) ; 



2/i(Q> M + - h(Q*){m N - m AQ ) + M&Xmff - mi) 7 „ 



:N(p) 



+ j.MQ') - f2(Q z )(m N + m Aq )| 7ai 4+2f 2 {Q')p^ 4 + |/ 2 (Q J ) + / 3 (<T) j (m w + m AtJ )<z„ 
+ |/2(g 2 ) + /3(Q 2 )U M ^-25i(Q 2 K75 + |5i(Q 2 )(^ + ^Aq) ~ 52(Q 2 )(m 2 v -m Ao )l 7M 7 5 



5i (Q 2 ) - 92(Q 2 ){m N - m AQ ) [■ 7a< #75 - 2g 2 (Q 2 )p f _ l #75 - ^ 52(Q 2 ) + 53(Q 2 ) f (m^v - m AQ )q^ 5 



92(Q 2 ) + g 3 (Q 2 )jq^ 4i5 J +•■• (9) 

In order to calculate the form factors /1, /2, /3, 31, 32 and 33, we will choose the independent structures p M , p^i 
q^i, Pnl5, Puils, and g M #y 5 from Eq. ©, respectively. 

For the theoretical side, to evaluate the correlation function in deep Euclidean region where (p + q) 2 -C 0, the 
explicit expression of the interpolating field of the Aq baryon is needed. Considering the quantum numbers, the most 
general form of interpolating current which can create the Aq from the vacuum is given as 

J Aq (aO = ^abc ^2(qf Cq\) l5 Q c + {3(qfC l5 q b 2 )Q c + (qfCQ b ) l5 q c 2 



+ P(qT C lb Q b )q c 2 + (Q aT Cql)j 5 ql + (3{Q a ' C 75 <72% \ , 



(10) 



where q\ and q 2 are the u and d quarks, respectively, a, 6, c are the color indices and C is the charge conjugation 
operator and (3 is an arbitrary parameter with (3 = — 1 corresponding to the Ioffe current. Using the transition current, 
J/I = <77/j(1 — 75 )Q and J AQ and contracting out all quark pairs by the help of the Wick's theorem, we achieve 



2(C) v< p(~f 5 ) p( 3 + {C) v0 (I)ps(j 5 )s<t> + {C)0<t>{l5\ 



+0 



2(C 7 5)^(/) P( 3 



(1 +75)7p S Q (-x)^{N(j,)\u a v (0)u b e (x)dl(0)\0), 



(11) 



where, Sq(x) is the heavy quark propagator which is given by [22 



5q(x) = (ac) - ig s 



2tt) 



ft + mg 



— 2 vx,G^ lv 



where 



S; 



47T 2 S 



(12) 



(13) 
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and Ki are the Bessel functions. Here, we neglect the terms proportional to the gluon field strength tensor since they 
can give contribution to four and five particle distribution functions and expected to be small [23L [24l [25[ . 

The matrix element (N(p) | e abc u^(0)u b g {x)d^(0) | 0) appearing in Eq. (fTTj) . which is the nucleon wave function, is 
represented as 0, H |M Hi H|| : 



4(0|e afcc <(«i^)4(a2a ; )^(a3^)l^(p)> 
= Sim N C a p{"ff ) N) 1 + S 2 m 2 N C a p(^ b N) 1 

+ Vim N { lb C) a0 N 1 + r 2 m 2 N ( l5 C) al 3^N) 7 + (Vi 



+ V 2 m N U>C) a0 (fa 5 N) 7 + V 3 TOjv(7MC)a/3 (7^75^)7 + V i m 2 N {j;C) a p{~f h N) 1 
+ V 5 m 2 N (~f^C) a p{ia^ x^N)^ + V 6 m^(^C*) a/3 (^75^) 7 + {Ai 

2 2 

mN{j>l§C) a p{fiN) 1 + -43"lAr(7Ai75C)a/3(7 /i ^ V )7 

+ A i m 2 N (ft^C) a pN 1 + Ar rn 2 N ('7 t _ l 'y b C) af )(i(j tl "x u N) 1 + Aem 3 N (^j 5 C) a( 3{^N) 7 

2 2 

+ (71 + ^LT^Wia^C)^^^ + T 2 m A r( 2; >^a AlI ,C)^(75A r )7 
+ %m N {(T liy C) a f3{cj^ l '^ 5 N) 1 +%m N {p u a lll ,C) a p{(T> ip x p ^ 5 N) 1 
+ %m 2 N {x v iaf lu C) a p(^^N) 1 +%rn 2 N {x^p u i(T^C) a f){^N) 1 
+ T 7 m 2 N (a^C) a p(a^^ 5 N) 7 +T 8 m%{x"a^C) a ^ p x p -f 5 N) 



7 > 



(14) 



where, the calligraphic objects which have not definite twists are functions of the scalar product px and the parameters 
ai, i = 1, 2, 3 and they are presented in terms of the nucleon distribution amplitudes (DA's) with definite and increasing 
twists. The scalar, pseudo-scalar, vector, axial vector and tensor DA's are explicitly shown in Tables HI HT1 IIII1 HVl and 
IVl respectively. 



Si = Si 



2pxS2 = Si — 5*2 



TABLE I: Relations between the calligraphic functions and nucleon scalar DA's. 



Pi =Pi 



2pxV 2 = Pi - P 2 



TABLE II: Relations between the calligraphic functions and nucleon pseudo-scalar DA's. 



Vi = 




2pxV 2 = Vi 


-V2-V3 


2V 3 = 


= v 3 


4pxV 4 = -2Vi 4 


■V 3 +V 4 + 2V 5 


ApxVs = 


V4-V3 


4(px) 2 V 6 = -Vi + V 2 


+ V 3 + V 4 + V 5 - Ve 



TABLE III: Relations between the calligraphic functions and nucleon vector DA's. 
The distribution amplitudes F(flipx)= Si, Pi, Vi, Ai, Tj can be written as: 

F{a t px) = J dx 1 dx 2 dx 3 6(x 1 + x 2 + x 3 - l)e- ipxSiXiai F(xi) . 
where, X4 with i = 1, 2, 3 corresponds to the longitudinal momentum fractions carried by the quarks. 



(15) 
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Ai = Ai 

2pxA 2 = -Aj + A 2 - A 3 

2A 3 = A 3 
ApxAi = -2Ai - A3 - Aa + 2A 5 

4pxA 5 = A3-A4 

4(px) 2 A 6 = Ai - A 2 + A 3 +A4- A 5 + A 6 



TABLE IV: Relations between the calligraphic functions and nucleon axial vector DA's. 



Ti — T\ 

2pxT 2 = Ti + T 2 - 2T 3 

2T 3 = T 7 

2pxTj = Ti - T 2 - 2T 7 

2pxT 5 = —Ti + T a + 2T 8 

4:(px) 2 T 6 = 2T 2 - 2T Z - 2Tj + 2T 5 + 2T 7 + 2T 8 

4pxT 7 = T 7 -T 8 

4(pxf% = -Ti + T 2 + T 5 - T 6 + 2T 7 + 2T 8 



TABLE V: Relations between the calligraphic functions and nucleon tensor DA's. 



In order to obtain the QCD or theoretical representation of the correlation function, the heavy quark propagator 
and nucleon distribution amplitudes are used in Eq. (jlip . Performing integral over x, equating the corresponding 
structures from both representations of the correlation function through the dispersion relations and applying Borel 
transformation with respect to (p + q) 2 to suppress the contribution of the higher states and continuum, one can 
obtain sum rules for the form factors fx, /b, fz, gx, 9i and 53. 

By means of the heavy quark effective theory (HQET), the number of independent form factors is reduced to two, 
Fx and F 2 . Hence, the transition matrix element can be parameterized in terms of these two form factors as (27l. [28j| : 

(N(p) I uTb I A Q (p + q)) = N(p)[Fx(Q 2 )+ ^F 2 (Q 2 )]Tu Aq ( P + q), 

(16) 

where, T is any Dirac matrices and f) = . One can immediately obtain the following relations among the form 

factors in HQET limit comparing the Eq. p^|) with the general definition of the form factors in Eq. © (see also 
[HHl) 

9l = fx = Fx + ^F 2 
F 2 

92 = h = .93 = h = ( 17 ) 

Considering the above relations, instead of giving the explicit expressions of the sum rules for the all form factors 
which are very lengthy, we will present only the expressions for fx and f 2 in the Appendix-A. However, we will give 
the extrapolation of all form factors in finite mass in terms of q 2 in the numerical analysis section. 

In the following, some remarks about how the HQET limit of the form factors satisfy the above relations are in 
order. In HQET, all the ratios, — , — , — , — , — , fy- and — should be equal to one. The deviation of those ratios 



from unity are presented in Tables rVTl and IVTll for A& — > piv and A c — > nlv, respectively. The bottom case and Lattice 
QCD input parameters satisfies the HQET relations exactly, while the maximum violation of this symmetry is related 
to the charm case and QCD input parameters. When we consider all relations, we see that the violations for charm 
case is larger than that of the bottom one. 

The explicit expressions of the sum rules for form factors reveals that to get the numerical values of the form factors, 
the expression for residue Aa q is needed. This residue has been calculated in [3l| using two-point QCD sum rules 
method: 

_ 2 

~X 2 AQ e~ mA Q / B = / e J ^p(s)ds + e^T, (18) 
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HQET 



QCD sum rulesinput parameters Lattice QCD input parameters 



fl 

91 







h, 

92 


20 


% 


h 

93 


20 


% 


h, 

93 


20 


% 


h. 

92 


20 


% 


f2 
f 3 







92 
93 








TABLE VI: Deviation of the ratio of the form factors from unity (violation of HQET symmetry relations) for Af, — > piu. 



HQET 

QCD sum rules input parameters Lattice QCD input parameters 



fl 

91 








h. 

92 


45 % 


40 % 


h, 

93 


45 % 


40 % 


Is. 

93 


45 % 


35 % 


h. 

92 


40 % 


35 % 


f2 

fa 








92 
93 









TABLE VII: Deviation of the ratio of the form factors from unity (violation of HQET symmetry relations) for A c —> nlv. 



with 



p(s) = (<rfd> + <^>)^-j|^[6(l + /J)Voo-(7+lWo2 

- 6(1 + /3)Vn] + (1 + 5(3)m Q (2^ w - ^ii - ^12 + 2^i)| 

+ ^dr^ 5 + + 5 0) W10 ~ 6^20 + 27/.30 - 4^4i + ^42 - 12/n(— )], 
ZU487T m,Q 



(19) 



^ < dd >K " u > 



22 9 



(!3 + "ffl + 77JT( 2 5 + 23,8) - (13 + 11/3) 



(20) 



where, so is continuum threshold, Ms is the Borel mass parameter and ip nm = 
functions. 



are some dimensionless 
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QCD sum rules [16] 


Lattice QCD [17, .18, 191 


In 


(5.0 ±0.5) x 10~ 3 GeV 2 


(3.234 ±0.063 ±0.086) x 10~ 3 GeV 2 


Ai 


-(2.7 ±0.9) x 10" 2 GeV 2 


(-3.557 ±0.065 ±0.136) x 10~ 2 GeV 2 


A 2 


(5.4 ± 1.9) x 10~ 2 GeV 2 


(7.002 ±0.128 ±0.268) x 10" 2 GeV 2 


vi 


0.23 ±0.03 


0.3015 ± 0.0032 ±0.0106 




0.38 ±0.15 


0.1013 ± 0.0081 ±0.0298 


ft 


0.40 ± 0.05 




fi 


0.07 ±0.05 




fi 


0.22 ±0.05 





TABLE VIII: The values of independent parameters entering to the nucleon DA's. The first errors in lattice values are statistical 
and the second errors represent the uncertainty due to the Chiral extrapolation and renormalization. 

III. NUMERICAL RESULTS 

The numerical analysis of the form factors and total decay rate for A h ( c ) — > p(n)£u transition are presented in 
this section. Some input parameters used in the analysis of the sum rules for the form factors are (uu)(l GeV) — 
(dd){l GeV) = -(0.243) 3 GeV 3 , m n = 0.939 GeV, m p = 0.938 GeV, m b = 4.7 GeV, m c = 1.23 GeV, m Ab = 
5.620 GeV, m Ac = 2.286 GeV and m§(l GeV) = (0.8 ±0.2) GeV 2 [Hj]. The main inputs which arc the nucleon DA's 
can be found in These DA's contain eight independent parameters /jv, Ai, A2, Vf , A", ft, ff and f£ . These 
parameters have been calculated in the light cone QCD sum rules [l6| and also most of these parameters have been 
computed in the framework of the lattice QCD [U El, [111 ■ For those parameters which have not calculated in lattice, 
the data from QCD input parameters will be used. These parameters are given in Table [Villi 

Three auxiliary parameters are encountered to the expression ofthe sum rules for form factors, continuum threshold 
so, Borel mass parameter Mj| and general parameter (3 entering to the most general form of the interpolating current 
for Aq baryon. A working region should be determined for these auxiliary and mathematical parameters such that 
the form factors as physical quantities should be independent of them. The continuum threshold, sq is not completely 
arbitrary and it is related to the energy of the exited states. From our results, we observed that the form factors are 
weakly dependent on s in the interval, {m^ Q + 0.5) 2 < s < {m\ Q + 0.7) 2 . To determine the working region for (3, 
we look at the variation of the form factors with respect to cosQ in the interval — 1 < cos6 < 1 which is corresponds 
to —00 < (3 < 00, where /3 = tanO. As a result, we attain a region at which the dependency is weak. The working 
region for (3 is obtained to be —0.75 < cosO < 0.25 for A& and —0.25 < cosd < 0.25 for A c . The Ioffe current which 
corresponds to cos8 = —0.71 is inside the working region for Af, but out of the region for A c . 

For further analysis, the upper and lower limits of M| should be determined. To do that, we apply two conditions: 
The first one, which gives the upper limit, is that the series of the light cone expansion with increasing twist should 
be convergent, and the second one, which determine the lower limit, is that the contribution of higher states and 
continuum to the correlation function should be enough small i.e., the contribution of the highest term with power 
l/Mg is less than, say, 20-25% of the highest power of Mg. In the present work, both conditions are satisfied 
in the region 15 GeV 2 < M 2 < 30 GeV 2 for A b and 4 GeV 2 < M 2 < 12 GeV 2 for A c , which we will use in 
numerical analysis. Taking into account the above requirements, we obtained that the form factors obey the following 
extrapolations in terms of q 2 : 

^ 2 )[.^ 2 )] = 7-Vl + 7T-V^' W 

y m 'jit 1 m 'fit 

The values of the parameters a, b and mfu for form factors and their HQET limit are given in Tables HXl [Xl IXII 
and IXIII related to the QCD sum rules and lattice QCD input parameters. Because of the working near the light 
cone, x 2 ~ and concerning the considered correlation function, the results are not reliable at low q 2 , hence to make 
the extension of our predictions to full physical region, we need to the above parameterization. From those Tables, 
we see that the pole of the form factors exist outside the physical region and the form factors are analytic in the 
whole physical interval. The values of form factors at q 2 ~ obtained from fit functions are shown in Tables IXIIII 
and IXIVl for A& — > "plv and A c — > nlv, respectively. A comparison of the existing predictions from other approaches 
is also presented for bottom case. The Table IXim depicts a good consistency on our result for /i(0) HQET limit 
obtained from lattice QCD input parameters with the prediction of [2l|, however the /i(0) HQET limit obtained from 
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QCD sum rules parameters is almost four times larger than that of [2l| prediction. On the other hand, the similar 
comparison of our result on form factor f 2 (0) at HQET and prediction of [2l[ shows that the value presented in [2l[ 
is almost two times greater than our result obtained from Lattice QCD input parameters and 1.5 times smaller than 
our result obtained from QCD input parameters. 





QCD sum rules 


Lattice QCD 




a 


b 


rrifit 


a 


b 


rrifit 


h 


0.025 


0.052 


4.91 


0.048 


0.016 


4.89 


h 


0.007 


-0.050 


4.92 


-0.003 


-0.006 


4.92 


h 


0.052 


-0.13 


4.99 


0.028 


-0.063 


4.96 


9i 


-0.059 


0.13 


5.29 


-0.17 


0.32 


5.32 


92 


0.011 


-0.050 


5.20 


0.019 


-0.040 


5.40 


93 


-0.009 


-0.017 


4.90 


-0.015 


0.012 


4.98 



TABLE IX: Parameters appearing in fit function of the original form factors for — > piv. 





QCD sum rules 


Lattice QCD 




a 


b 


rrifit 


a 


b 


rrifit 


fl 


-0.034 


0.20 


1.59 


-0.14 


0.64 


1.55 


h 


-0.015 


-0.77 


1.57 


0.018 


-0.32 


1.60 


h 


-0.062 


-1.23 


1.48 


0.12 


-1.09 


1.56 


9i 


-0.015 


0.54 


1.53 


-0.20 


0.71 


1.59 


92 


-0.11 


-0.20 


1.52 


-0.034 


-0.14 


1.65 


93 


-0.088 


0.085 


1.48 


0.009 


-0.41 


1.50 



TABLE X: Parameters appearing in fit function of the original form factors for A c — » nlv. 





QCD sum rules 


Lattice QCD 




a 


b 


rrifit 


a 


b 


rrifit 


fl 


0.041 


0.040 


4.82 


0.0042 


0.016 


4.92 


h 


0.033 


-0.097 


4.83 


0.013 


-0.030 


5.92 


h 


0.060 


-0.14 


4.90 


0.016 


-0.040 


4.94 


9i 


-0.0012 


0.096 


5.10 


-0.0022 


0.029 


5.30 


92 


-0.0094 


-0.018 


5.36 


0.0017 


-0.0043 


5.36 


93 


-0.040 


0.025 


4.95 


-0.018 


0.015 


4.98 



TABLE XI: Parameters appearing in the fit function of the form factors at HQET limit for At — > piv. 

In the next step, we calculate the total decay rate of Aq — ► Ntv transition in the whole physical region, i.e., 
mf < q 2 < (m^ Q — m^) 2 . The decay width for such transition is given by the following expression [33l. |34| 

A 2 

T(A Q - Nln) = 384 ^ 3 \V qQ \ 2 J dq 2 (1 - mf/q 2 ) 2 ^/ (Z 2 - q 2 )(A 2 - q 2 ) N(q 2 ) (22) 

Q 9 
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QCD sum rules 


Lattice QCD 




a 


b 


rrifit 


a 


b 


rrifit 


fl 


-0.066 


1.14 


1.51 


-0.039 


0.37 


1.55 


h 


0.046 


-1.14 


1.53 


0.047 


-0.63 


1.48 


h 


0.071 


-1.33 


1.50 


0.039 


-0.52 


1.53 


9i 


-0.10 


1.21 


1.57 


-0.039 


0.39 


1.55 


92 


-0.070 


-0.11 


1.56 


-0.027 


-0.046 


1.60 


93 


-0.076 


-0.91 


1.54 


-0.034 


-0.032 


1.54 



TABLE XII: Parameters appearing in the fit function of the form factors at HQET limit for A c 



nlv. 





Original 


HQET 




QCD sum rules 


Lattice QCD 


QCD sum rules 


Lattice QCD 


[21] 


/i(0) 


0.077 


0.064 


0.081 


0.021 


023+ ' 006 

u - uio -0 .005 


/a(0) 


-0.044 


-0.013 


-0.064 


-0.018 


-0 39+ 006 

u - UJ3 -0 .009 


MO) 


-0.079 


-0.036 








9i (0) 


0.073 


0.15 








52(0) 


-0.039 


-0.021 








03(0) 


-0.026 


-0.0035 









where 



TABLE XIII: The values of the form factors at q 2 = for A 6 



Ntf) = W)(AW -mi) + 2£^A 2 (l + 2m 2 /«^) - + 2«^)(2 (? 2 +mf)) 

+ F 2 (q 2 )(A 2 - g 2 )(2£ 2 + q 2 )(2q 2 + m 2 )/m| b + 3F 3 V)m 2 (£ 2 - <rV/m 2 b 
+ 6F 1 (q 2 )F 2 (q 2 )(A 2 - q 2 )(2q 2 + m 2 )£/m Efc - m{q 2 )F^)m 2 {Y? - q 2 )A/m^ b 
+ G 2 (q 2 )(Y?(Aq 2 - m 2 ) + 2S 2 A 2 (1 + 2m 2 /g 2 ) - (A 2 + 2q 2 )(2q 2 + m 2 )) 
+ G 2 (g 2 )(S 2 - 9 2 )(2A 2 + q 2 )(2q 2 + m 2 )/m| b + 3G 2 ( 9 2 )m 2 (A 2 - q 2 )q 2 /ml b 
- 6G 1 ( (7 2 )G 2 (g 2 )(E 2 - q 2 )(2q 2 + m 2 )A/m Sb + 6G 1 ( g 2 )G 3 ( g 2 )m 2 (A 2 - <? 2 )£/m Sb 



(23) 



Here, F^q 2 ) = h{q 2 ), F 2 (q 2 ) = m AQ f 2 (q 2 ), F 3 (q 2 ) = m AQ f 3 (q 2 ), G 1 (q 2 ) = 9l {q 2 ), G 2 (q 2 ) = m AQ g 2 (q 2 ), G 3 (q 2 ) = 
m AQ g 3 {q 2 ), S = m\ Q + mn and A = m\ Q — mjv- Gf — 1.17 x 10 -5 GeV -2 is the Fermi coupling constant, and mi 
is the leptonic (electron, muon or tau) mass. For the corresponding CKM matrix element V u b = (4.31 ± 0.30) 10~ 3 
and Vd = (0.230 ± 0.011) are used [35 1. 

Our final results for total decay rates are given in Table IXVI As it can be seen from this Table, our results for e 
and fi and A& cases are consistent for two sets of input parameters when the original form factors are used especially, 
when we consider the uncertainties. However, QCD input parameters result is 1.5 times greater than that of the 
lattice input parameter for the decay rates of r and bottom case. If we consider A c , QCD sum rules input parameters 
gives the result 2 times greater than the Lattice QCD input parameters. On the other hand, when we consider the 
uncertainties, results obtained using both sets of input parameters and original form factors coincide for all leptons. 
At HQET limit and QCD sum rules input parameters, our predictions for the decay rates are in the same order of 
magnitude with the original form factors and two sets for all leptons and both charm and bottom cases. In contrast, 
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Original 
QCD sum rules Lattice QCD 


HQET 

QCD sum rules Lattice QCD 


/i(0) 


0.17 


0.50 


1.078 


0.33 


/a(0) 


-0.78 


-0.31 


-1.09 


-0.58 


h(P) 


1.29 


0.98 






9i (0) 


0.52 


0.51 






92(0) 


-0.31 


-0.18 






33(0) 


-0.0032 


-0.31 







TABLE XIV: The values of the form factors at q 2 = for A c -> nlv. 





A fa ► PM^fi 


A fa ► pei , e 


A b P TU T 


A c «■ npw^ 


A c • nev e 


A c ► nrv T 


For QCD sum rules inputs 


(3.07 ± 1.05) X 10 -16 


(3.065 ± 1.05) x 10 -15 


(3.82 ± 1.35) X 10~ 15 


(2.89 ± 0.95) X 10~ 13 


(2.86 ± 0.95) X 10~ 13 




For lattice QCD inputs 


(2.87 ± 0.95) X 10 — 15 


(2.87 ± 0.95) x 10 -15 


(2.55 ± 0.85) X 10~ 15 


(1.35 ± 0.45) X 10 -13 


(1.33 ± 0.43) X 10 -13 




HQET limit for QCD sum rules inputs 


(5.84 ±1.81) X 10 — 15 


(5.83 ± 1.81) X 10 — 15 


(7.90 ± 2.45) X 10~ 15 


(5.08 ± 1.65) X 10~ 13 


(5.01 ± 1.60) X 10~ 13 




HQET limit for lattice QCD inputs 


(4.70 ± 1.60) x 10 — 17 


(4.60 ± 1.55) x 10~ 17 


(2.36 ± 0.85) X 10~ 16 


(8.75 ± 2.85) X 10~ 14 


(8.74 ± 2.83) X 10~ 14 





TABLE XV: Values of the F(Aq — > Nlv) in GeV for different leptons and two sets of input parameters obtained from QCD 
sum rules and lattice QCD and also their HQET limit. 

the results at HQET limit and lattice parameters are two orders of magnitude less than HQET limit and sum rules 
inputs as well as original form factors for bottom and e and \x cases. For r and bottom, and e and fj, and charm 
cases, this difference is approximately one order of mag nitude. We also compare our results on decay rates in units of 
\VqQ\ 2 s' 1 with the predictions of references [1, OH [LIT E3, HE HH in Table IXVll From this Table, it is clear that our 
results for bottom case, lattice parameters and HQET limit are in the same order of magnitude with the predictions 
of [HI, and HQET-[20l|. For all other cases the difference between our results with the existing predictions of the 
other approaches presented in Table IXVII is one-two order of magnitudes. In Table I XVII HOSR refers to harmonic 
oscillator semi relativistic and HONR stands for harmonic oscillator non relativistic constituent quark models. 

To summarize, using the most general form of the interpolating currents of Aq and nucleon DA's with two sets 
of input parameters, namely QCD sum rules and lattice QCD inputs, the transition form factors of the semileptonic 
Aq — > Nlv have been calculated in the framework of the light cone QCD sum rules in full theory and HQET. The 
lattice input parameters satisfy the HQET relations exactly for bottom case, while the maximum violation is for 
charm case and QCD input parameters. The results of the form factors at HQET and q 2 = have been compared 
with the existing predictions of the other approaches. These transition form factors have been used to estimate 
the corresponding tree level semileptonic decay rates both in full theory and HQET limit. A comparison of the 
obtained results and the existing predictions of the other approaches which all are at HQET limit, was also presented. 
The best consistency between our results and those predictions is related to the bottom case and lattice QCD input 
parameters at HQET. Our results can be checked in experiments hold in future such as LHC. Comparison between the 
experimental data and our results could give essential information about the nature of the Aq, Nucleon distribution 
amplitudes as well as determination of the CKM matrix elements, V u b and V c d- 
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A b PM^JLt 


A b ► pei/ e 


A b ► p-rv-r 


A c ► nf.iu^ 


A c > ne 




Ac ► n-rv-r 


For QCD sum rules 


(2.5 ± 0.85) X 10 14 


(2.5 ± 0.85) X 10 14 


(3.12 ± 1.05) X 10 14 


(8.3 ± 2.85) X 10 12 


(8.21 ± 2.80) 


< 10 12 




For lattice QCD 


(2.35 ± 0.85) X 10 14 


(2.35 ± 0.85) x 10 14 


(2.08 ± 0.70) X 10 14 


(3.88 ± 1.25) X 10 12 


(3.82 ± 1.20) 


< 10 12 




HQET limit for QCD sum rules 


(4.78 ± 1.75) X 10 14 


(4.77 ±1.75) x 10 14 


(6.46 ± 2.15) X 10 14 


(1.46 ± 0.55) x 10 13 


(1.44 ± 0.55) 


< 1013 




HQET limit for lattice QCD 


(3.84 ± 1.25) x 10 12 


(3.76 ± 1.20) X 10 12 


(1.93 ± 0.70) X 10 12 


(2.51 ± 0.85) X 10 12 


(2.51 ± 0.85) 


< 10 12 




[10J 


2.05 x 10 13 






2.58 X 10 13 






6.48 X 10 12 




QCD sum rulcsf20j 


3.65 X 10 13 




HQET[20] 


5.62 x 10 12 




[121 


4.55 X 10 12 (HONR) 


4.01 X 10 12 (HONR) 


1.02 x 10 10 (HONR) 








7.55 X 10 12 (HOSR) 


6.55 X 10 12 (HOSR) 


1.35 X 10 10 (HOSR) 







TABLE XVI: Values of the total decay rate (in IVqQj 2 s 1 ) of the Aq — > Nlv transition for different leptons and two sets of 
input parameters obtained from QCD sum rules and lattice QCD and also their HQET limit compared to the 0, Hol . [Til . Il2l 
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Appendix A 



In this Appendix, the explicit expressions for the form factors /i and fa are given: 
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+M B t ^/3W-i2 2 , + 16 6 , + 18,+20 19 ,+22 24 ,-24 48 ,+6 2 (x») + toW-12 4 ,-18,-20i 3 ,-22 24 ,-6 4 ( x i) 

+to/3W + i2 4 ,-i8 3 ,-2033 > -22 12> +6 4 (^)) J + jg 4 (-l + *o)(l + 2tg) + g 2 M 2 t (-1 + to - 2i§ - 3tg + 3tg) 



^(2Q 2 + 3M|t )s(s ,Q 2 



(1-/3) 6Wi (a:i)(to - x 2 ) + 2x 2 H 16 (x t ) - 4x 2 H 24 (x l )(5 + 4/3) 



+4H 22 (x l )(2 + /3) 



+tg(-2 + t )(2g 2 + 3M 2 B t )s{s ,Q 2 



Q 2 Ml(-1 + t )(2io - 2t 2 - tg + 9t 4 ) + Q 2 (-2 + 4t - t 2 - 4tg + 6t 4 ) 



M 2 t 



+M 2 Q 2 (-2t + 3t 2 -t 3 - 7t 4 + m 5 Q ) + 4{-l + t )(2Q 2 + 3t M 2 )s(s , Q 2 ) + Q 2 t 2 



+ 



— Q 2 tuH + iQ-2i 2 { x i) 



-t 2 ,M B W+i2 2 , + 14, + 15, + 15,-204,+21 4 -24 2 ,-6 2 ,+8,+9(xi) + to^io (Xj) (Q 2 + to-M s ) — X 2 (Q 2 + to M 2 J )W + i ,-l6,+24 2 { x i) 

| - MlQ 2 (-l + /3)W +122 ,_ 182: _ 62 (x l ) - M 2 Q 2 W 20 (x,)(2 - 22/3) + M 2 Q%(-1 + P)H +12l0 ,-is,+6 e ( x i) 

— MgQ 2 t 'H2o(Xi)(23 + 85/3) + g 4 to/37Y+io 6 , + 16 2 ,+22 8 ,-24 16 (xi) + Q 4 £ 2 ,^-10 6 ,-16 2 ,+22i 6 ,-24 2 o(Xi) 
+M%Q 2 t 2 t H- We , + 12-8, -14, +155,-16 2 ,+18 5 ,-2 8 ,+20 2 i,+21 20 , +23 36 , -24 20 ,+4 4 ,+6 8 ,+83,-9 7 (Xj) 
+M|(5 2 t 2 /3H + io 6 ,-l 28 , + 14,-155, + 16 2 , + 187,-4 2 , +20 6 i, +21ie, +23 5 6, -24 16 ,+4 8 ,-6 8 , -8 3 , +9 7 (xi) 
+M|Q 2 tQ , H_10 6 , + 124,-16 2 ,-186,-20 6 ,+223 2 ,-24 2 0,-64 (Xi) + Q ^oH+lo 6 , + 16 2 ,-22 24 ,+24 20 ( x i ) 
+^iQ 2 io/3W+10 6 , + 124, + 16 2 , + 18 2 ,+2038,+22i 6 ,-24i6,+6 4 (Xi) + Q 4 ^o/3H-lo 8 ,-16 2 ,-22i 2 ,+24 16 ( x i) 

10 6 ,-12i 2 , + 16 2 , + 18 2 ,-20 22 ,-22 24 ,+24 2 o,-6 4 (Xi) + M B Q t (3H- 10 6 , + 12i 2 ,- 16 2 ,- 1810 ,-20no ,-22i 2 ,+24 16 ,+64 { x i ) 



+s(so,<9 2 



Q 2 to^ -10 6 ,-16 2 ,+22i 6 ,-24 20 , + 16 2 ,+22 8 ,-24i 6 (x i )-4Q^W 22 (x i )(2 + /3) 



+M B t W_io 6 ,+12 4 ,-16 2 ,-18 6 ,-20 6 ,+22i 6 , -2420,-64(^1) + M B i o/^+106,-124, + 16 2 , + 18 2 ,+20 38 ,+22 8 ,-24 16 ,+6 4 ( x 



+M B i 4 7Y_i2 8 ,_i8 2 ,-20 26 ,-22 8 ,+6 8 (xi) + M s to/3W+i2 8 ,-18 6 ,-20 66 ,-22 4 ,+6 8 (Xi) 



2+4, 



+ 2t x 2 



2Q 2 W 22 (x,)(2 + /3) 



M 2 t (-4 + 3f ) + Q 2 (-2 + 3t ) - 2Q 2 H 2 4(xi)(5 + 4/3) Q 2 (-l + t ) + M 2 (-1 - t + t 2 ) 



+2 H 24 (x t )(5 + 4/3) + H 22 (x t )(2 + /3)(-2 + t ) (Q 2 + M 2 t )s( Sa ,Q 2 ) + W +10 „i 6 (a;i)(-l + /3) 



Q 4 (-l + 1 ) + Q 2 M 2 (-1 - 1 + *o) - (Q 2 + M 2 t )s(s , Q 2 ) 



m 2 N m b Q 2 



M 2 



(-1+/3) Q 2 t M 2 



(-l + to)to 



W+i2 2 ,+i4,+i5,-20 4 ,+2i 4 ,-6 2 ,+8,+9(x 4 ) + 2Q%H 22 (xi) [Q'(2 - 3t ) + 6M B t (l - to) J - 2W 22 (a; 
(g 2 + M 2 t )s(s , Q 2 ) - x 2 Q 4 (-l + 3t ) + x 2 (-l + 1 ) (V 2 t s( So , Q 2 ) + Q 2 (3M 2 t + s(s , Q 2 )) 
+m N Q% |g 2 ^2W 6 (xi)(-l + /3)(1 + to) + Wig (a*) (-3 + /3 + t - 5t /3) - W 2 o(xi)(3 - 19/3 + llt + 55t /3)^) 
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2H 6 (x l )(-l + /3)(1 + 2t ) + W 20 (^)(3 - 19/3 + 13t + 33t /3) + Hi 8 (x,)(3 + t - /3 + 3t /3) s(s , Q 2 ) 



+2Wi 2 (a;i)(-l + 0) Q 2 (-l + 3t ) + (-1 + 2t )s(s 0) Q 2 ) 



+ 



:(*0 - ^2) S - TOAT^O 



(g 2 +m2 r t2)2 2V 6^ 

(-1 + /?) ( m|r(-l + to) + m N m b t Q + Q 2 t (-1 + t ) - t s(s 0> Q 2 ) ) + Wio(xj)(-l + (3) [ 3m 2 N (-l + t ) - m N m b t 



+3i (0 2 (-l + to) - s{s , Q 2 ))^j - 2H 2i (x l ) (m%(5 + 4/5)(-l + t ) + m 2 N m b (-l + /3)t + (5 + 4/3)t (V(-l + to) 



-s(s ,Q 2 ) 



m| r m 6 (-l + /?)(-2 + t )(-l + t ) + 2m%(2 + /3){2t - 4t 2 + tg) - 2to w (2 + /3) 



+ H 22 (Xi) 

I'o [ Q 2 l -2 + 3t ) + (-2 + toHs-o.C^) ) +'»/,(.- I + <)/.. I ^ J i- i - 3/,,) -- (-I - /o .)«(*„. 

TOAf 



2mNH22{xi){to - x 2 ) 



m z N m b (-l + /?)(-2 + to)(-l + t ) + 2m Ar (2 + /3)t 



(g 2 + m 2 ^t2)4 v / 6M| to 

(2 - 4t + t 2 ,) - m 6 (-l + /3)to (V + 2t M| - 3g 2 t + s{s , Q 2 ) + t s(s a ,Q 2 )^j + 2m N {2 + /3)t 2 (q 2 (2 - 3t ) 



+2t Mj, + 2s( So , g 2 ) - t s(s , g 2 ) 



+ to 



- 2m 2 N MlHis.+20^+e(xi) + m 2 N M B (3Hi8 2 , +2 o 22 ,+6 2 ( x i) 

— m 4 N t H + io 6 , + W 2 ,+20 2 4( x i) + m Ar-^BtoW + i8,-2023,-66 ( x i) + m NTn b MgtoH + i^ + 15,-20i,+21 4 ,-6 2 ,+S,+9(xi) 

— MgQ 2 t a H+i$ 3 ,-20 3 ,-6 2 ( x i) + fn%PtoT^+W 6 ,+16 2 ,-24 16 (Xi) — m 2 N Mgt PH-is 7 ^20 ss ,+6 6 ( x i) 

-m N m b Mgt l3H + i4 i+ i5-2o 4 ,+2i 4 -6 2 ,+8,+9(xi) + M B Q 2 t /3H + i8,+20 19 ,+6 2 {xi) + m 2 N tl(m 2 N - Q 2 ) 
W+io 6 ,+i62,+24 2 o( x i) — m3 N m bt 2 H+\o 2 ^i& 2 ,+2A l {xi) + m 2 N (3tQ{Q 2 — m A r)W + 10 6 ,+16 2 ,-24i 6 (^i) 

— m N M B t H+10 e ,+14,-15 5 , + 16 2 ,-18 5 , +82,-2021, +2l2o,-23 3 6,+24 2 o,-44, -6 8 ,-8 3 , +9 7 ( x i) + m 3 N m b (3tlH + io 2 ^16 2 ,+24 4 (xi) 
+M 2 B Q 2 t 2 ) Ti + ia-20 11 -& 2 {xi) + TO 2 v M|t§/37Y + i0 6 , + 14,-155, + 16 2 , + 18 7 ,-2 4 ,+20 6 i,+21i 6 ,+2356,-24i 6 ,+4 8 ,+6 8 ,-8 6 ,+9 7 (a;i) 
-M B Q f3t H + ls 5 ,+20 5 5-62( x i)+ m NQ t H + 10 6 , + 162,+24 2 o (Xi) - m 2 Vg 2 /3ti5H + io 6 ,+i62,-24i 6 (a; 4 ) 

+.s(s , g 2 ) — M B to W + 18 3 ,+20 3 , +62(^1) + ^i/?toW+18,+20i9,+6 2 ( a; i) _ m Arto^+10 6 , + 162 ,+24 20 ( x i ) 
-MstlTi+is, +2O13, +64 ( x i) + m Art 2 /3H+io 6 , + 16 2 ,-24i 6 (a;i) - M B tl(3H + 18 3l +2033,-64(^0 

-2M!Wi 2 (zi)(-l + /J) ^TOArTObto + m 2 N (-l + t )(-l + 4t ) + t g 2 - 3Q 2 tg + t s(s , g 2 ) - 2t 2 s(s , g 2 ) 
+2x 2 m%n w (xi)(-l + f3) (m 2 N (-l + t ) + m N m b t - t M| - t a Q 2 + Q 2 t 2 - t s(s , Q 2 ) 

+2x 2 m 2 N H w (x l )(-l + [3) (^m 2 N {-\ + t ) - m N m b t - 3t (M 2 + g 2 (l - t )) - 3t s(s , g 2 )J + 4x 2 m 2 N H 24 (xi) 
m 2 N (5 + 4/3)(-l + t ) - m N m b (-l + f3)t + (5 + 4/3)t (V 2 + g 2 (l - f ) + s(s , Q 2 )^j ) } } 
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and 

h{Q 2 ) 



1 /-1-X2 



]__ e ml b /M 2 B 



dx 2 

to J a 



i—x 2 



2V6x 2 



(-1 + f3)H+u >+5 (xi) + H+n(xi)(-l + 6/3) 



r 1 fl—X2 PX2 „ „ 'I 

+ dx 2 dx x \ dt x e- s ^W M » 

Jto Jo Jto \M%tl 2 2s/Q 



- m N m b (-l + (3)x 2 H+i ,-i6,+24 2 (^i) + 2W 22 (xi) 



m N m b (-l + /3)x 2 + (Q 2 + s(t u Q 2 ))(-l + /3)x 2 + m 2 ^! - /3 + /3a; 2 + hx 2 ) 



+ 



M|tl2v/6 . 



m N m b (-l + (3) 



W +10 ,-i6,+24 2 (a; i ) - 2H 22 (x,) ( m iV m 6 (-l + /?) + (Q 2 + s(i i5 Q 2 ))(-l + /3) - ( 3Q 2 + s(ti, Q 2 )(2 + 0) ) x 2 

™2 



+m^(5 + /3 + 2a; 2 +/3x 2 ) 



JV 



M|*14v/6 



14, + 155, + 186,-82,+20io,+2l2o,-224,+3623,+4 4 ,+64,+83,-97 



0»i) 



-AHi 2 {xi)( — 1 + /?) + /3W + i4 ; _i5 5j+ i8 6 ,-24,+2034,+21i6,-22 8 ,+4 s ,-64,-83,+9 7 ( a; i) — 8x 2 (2 + f3)H 22 (Xi) 



+ 



M4n 3 V6 



^2(-l + /3)W 22 (x 4 ) + 



M2il22V6 



m 6 (-l + /?)W_i2,+i8,+20,6(a;i) + 2m N x 2 (l + 2(3)H 22 (xi) 



2 '" 2a ' (2 + 



+ dx 2 d Xl e- So/M z { (t - x 2 ) 

Jto Jo j 



171 N^O 



N 



m N m b t (-l + f3) 



(Q 2 + m 2 N t 2 fV6 (Q 2 + m 2 N t 2 ) 2 2V& 
H +w ,-i6,+24,2( x i) + 2H 22 (x l ) ^ - m N m b t (-l + (3) - Q 2 t (-1 + (3 + 3to) + m 2 , ^1 - (3 - i„(5 + /?) + (2 + /3)i 2 



to -l + (3+(2 + f3)t s(s ,Q 2 ) 



+ 



m N 



4:m N H 22 (xi) (t -x 2 )\ - m N m b (-l + (3)t 



(Q 2 + m 2 N t 2 )M 2 t 4V6 

+m^(l -0-5to- 0k + 2t 2 Q + (3t 2 Q ) + to (q 2 - Q 2 (/3 + 3t ) + M§(— 1 - 2/3 + 2/3t + 4t ) + *(«o, Q 2 ) 
-(/3 + 2t +/3t )s(so,<3 2 ) ) ) + *o \m b M] 3 H- 2 02 -62,+2o 2 ,+6 2 (xi)+m 2 N m b to7i-io2,+i62 -24 4 (xi) 



+M 7 2 jm6toW-14, + 155,-2 8 ,+20io,+2l20,+2336,+4 4 ,+6 4 ,+83,-97( x i) + m N m b (3to'H+102 ,-16 2 ,+24 4 ( x i) 

+M%m N [3toH + iA ,-i5 5 ,-2 4 ,+2034,+2ii6,+2356,+48 ,-6 4 -83,+9 7 (zi) _ 2M%H\ 2 (x^ (- 1 + /3) (m b + 2m N to) 
+2M B H ls (x l ) (m b (-l + f3) + 'Am N to(l + 0)j - 2m 2 N m b x 2 (-l + (3)H +W -i6,+24 4 ( x i) 

where 



(A2) 



W(a;i) = H(xi,x 2 , 1 - £1 - z 2 ), 
S (y,Q 2 ) = (l-y )m 2 N + ^—^Q 2 +™l, 

v y 

and to = to(so, Q 2 ) is the solution of the equation s(to, Q 2 ) — s , and is given as 
,2 _ n 2 



toOo 



2 _ m 2 N -Q 2 + y/4m 2 N ( m j + Q 2 ) + (m 
1 ~ 2 



so) 2 + «o 



2m 2 N 



(A3) 



(A4) 



In calculations, the following short hand notations for the functions H±i at ±j bt ... — ±aHi ± bHj... have been used, 
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and Hi are given in terms of the DA's as follows: 

Hi = Si H.2 = 5*1,-2 

H.3 = P\ H4 ~ Pl,-2 

H5 = Vi Tie = ^1,-2,-3 

n 7 = V 3 H & = -2^i,_b + ^3,4 

= ^4,-3 W10 = —14,-2,-3,-4,-5,6 

Wn = Ai H12 = — j4i,_2,3 

W13 = ^3 W14 = — 2Ai,_5 — A 3 ^4 

W15 = ^3,_4 Wi6 = ^1,-2,3,4,-5,6 

W17 = Ti His = 7l,2 — 2T 3 

Til9 = ^7 ^20 = ^1,-2 — 2T7 

W21 = — Ti,_5 + 2T§ W 2 2 = 72,_3,_4,5,7,8 

^23 = ^7,-8 ^24 = —^l, -2,-5,6 + 2T 7 , 8 , (^-5) 

where for each DA's, A±j,-y,... = ±Xj ± A.,-... have also been used. 



